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ABSTRACT 

We present 3D resonant radiative transfer simulations of the spatial and spectral dif- 
fusion of the Lya radiation from a central source in the host galaxies of high column 
density absorption systems at z ~ 3. The radiative transfer simulations are based on a 
suite of cosmological galaxy formation simulations which reproduce a wide range of 
observed properties of damped Lya absorption systems. The Lya emission is predicted 
to be spatially extended up to several arcsec, and the spectral width of the Lya emis- 
sion is broadened to several hundred (in some case more than thousand) km s -1 . The 
distribution and the dynamical state of the gas in the simulated galaxies is complex, 
the latter with significant contributions from rotation and both in- and out-flows. The 
emerging Lya radiation extends to gas with column densities of Nm ~ 10 18 cm -2 and 
its spectral shape varies strongly with viewing angle. The strong dependence on the 
central Hi column density and the Hi velocity field suggests that the Lya emission will 
also vary strongly with time on timescales of a few dynamical times of the central re- 
gion. Such variations with time should be especially pronounced at times where the host 
galaxy undergoes a major merger and/or starburst. Depending on the pre-dominance of 
in- or out-flow along a given sightline and the central column density, the spectra show 
prominent blue peaks, red peaks or double-peaked profiles. Both spatial distribution and 
spectral shape are very sensitive to details of the galactic wind implementation. Stronger 
galactic winds result in more spatially extended Lya emission and - somewhat counter- 
intuitively - a narrower spectral distribution. 
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1 INTRODUCTION 

Searches for Lya emission that are performed down to very 
low flux levels can detect very faint objects where the Lya 
emission is little affected by dust. Lya emission has there- 
fore dev eloped into an important tool to study high-redshift 
galaxies JCowie & H u 1998; St eidel et al j200d:lHu et al .1120021 ; 
Kodaira et all 120031: iRhoads et al 120041: iTaniguchi et"ai1 | 2005| : 
Kashikawa et al l 120061: live et al. 120061 : IStanwav et all 120071 : 
Ota et al.ll2008l : buchi et al.ll201Qh . At high redshift, or for in- 
trinsically faint objects, it is also often the only means of es- 
tablishing a spectroscopic redshift. At low flux levels, Lya 
emission due to star formation in galaxies can be spatially 
very extended due to the resonant scattering of Lya pho- 
tons — even if Lya cooling is not important and the corre- 
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spond i ng continuum emission is very comp act ( Havashino et al.l 
2004; iDiikstra. Haiman, & Soaansl 120061 : iRauch et all 12008: 
Finkelstein et al.ll2010h . 



Deep Lya observations of the high-redshift Universe thus 
have tremendous potential to probe the surrounding neutral gas 
in protogalaxies, and thereb y shed light on ga laxy form ation . 
The spectroscopic survey of IRauch et al.1 J2008l hereafter [R08h 
for low surface brightness Lya emitters based on a 92 hour 
long exposure with the ESO VLT FORS2 instrument has pushed 
Lya surveys to new sensitivity limits and yielded a sample of 
27 spatially extended (several arcsec) low-surface Lya emit- 
ters with fluxes of 3 x 10" 18 ergs -1 cm -1 and a space density of 
~ 3 x 10" 2 /^ Mpc -3 , similar to that of the faintest dwarf galax- 
ies in the loca l Uni verse. Based on the inferred incidence rate 
for absorption jR08l argued that their sample should overlap sig- 
nificantly with the elusive host population of Damped Lya Ab- 
sorption Systems (DLAs). 



© 2010 RAS 



2 Barnes, Haehnelt, Tescari and Viel 



.Barnes & Haehneltl (120091 . I2O10L hereafter iBHOl 
BHICh. building on the su cc essful model for DLAs of 



Haeh nelt. Steinmetz. & Rauchl dl998l . l2000h . presented a 



simple model that simultaneously accounts for the kinematic 
properties, column density distribution and incidence rate of 
DLAs and the luminosity function and the size distribution 
of the !r08| emitters in the context of the ACDM model for 
structure formation. The modelling of lBHld assumed spherical 
symmetry and simple power law profiles for gas density and 
peculiar velocity. Lya radiative transfer is, however, sensitive 
to the details of the spatial distribution and velocity of neutral 
hydrogen — inhomogeneity can provide channels for rapid 
spatial diffusion, while bulk velocities have a similar effect in 
frequency space. 

The complicated physics of galaxy formation — gravity, 
hydrodynamics, radiative cooling, star formation, supernovae, 
galactic winds, ionising radiation, metal and dust production — 
has led to numerical simulations being the method of choice 
for investigating the physical properties of the gas from which 
galaxies form. These simulations provide us with realistic 3D 
density and velocity fields that take into account most of the 
relevant physics. 

He re we w ill apply the Lya radiative transfer code devel- 
oped in IbHIOL adapted for the post-processing of 3D gas dis- 
tributions, to the results of simulation of DLA host galaxies in 
a cosmological context by iTescari et" ID {2009). We study the 
effect of the spatial distribution and the kinematic state of the 
neutral hydrogen on the low surface brightness spatially ex- 
tended Lya emission, which should become an important probe 
of kinematics of the gas in DLA host galaxies. 

In this first paper we will focus on a detailed investiga- 
tion of a small number of such simulated systems. In Section 
[2 we describe the (minor) changes to our Lyor radiative trans- 
fer code needed to follow photons through a density distribution 
on a 3D grid. In Sectio n IXTI we describe the set of simulations 
of ITescari et al ] d2009h that we have used, and in Section 13.21 
we describe the three haloes that we have extracted from the 
simulations for further study. Section HI discusses the haloes in 
detail, looking particularly at the effect of galactic winds on the 
gas properties and Lya emission. We discuss our results in Sec- 
tion[5] including their relevance to observed populations of high 
redshift Lya emitters, including LBG's, before presenting our 
conclusions in Section[6] 



2 THE 3D Lya RADIATIVE TRANSFER CODE 



been studied 


j a number of authors dCantaluDO et al.l 


2005; 


Tasitsiomill2006 


uLaursen. Razoumov. & Sommer-Larsenl 


2009; 


Zhens et al.ll200 


9tlKollmeier et al.ll2010l:lFaucher-Gi2uere et al. 


2010b. in the context of. for example, fluorescent emission from 



the IGM, bright Lya emitters at z « 8, young Lyman break 
galaxies and Lya emission linked to radiative cooling. 

We have previously described our imp lementation of a ID 
Lya radiative transfer algorithm in iBHIOl . In short, Lya pho- 
tons are created according to the emissivity profile, and then 
propagated through the gas. Each time the photon is scattered, 
a new frequency is calculated by assuming that the photon's 
frequency in the scattering atom's frame is unchanged. Thus, 
the code tracks the random walk in both real and frequency 
space before the photon escapes and is observed. The new code 
used here follows photons through a regular 3D grid. The ac- 



celeration scheme, by which photons which have frequencies 
close to line centre are scattered into the wings, is adapted 
to the conditions in the current cell. We have tested and sub- 
sequently adopted a prescription simil a r to t hat described in 
lLaursen. Razoumov. & Sommer-Larsenl (12 009). Further details 
on the code can be found in lBarnesI (120101) . 

The emergent spectrum will depend on the angle 
from which the system is viewed. Rather than "wait- 
ing" for enough photons to emerge in a given direc- 
tion, we implement the "peeli n g-off" algo rithm, described in 
lYusef-Zadeh. Morris. & White] dl984l) and IWood & Revnoldsl 
(Il999l) . At each scattering, the probability of escape in the direc- 
tion of the observer is calculated and a suitable weight added to 
the corresponding 2D pixel. The photons that eventually escape 
the system are used to calculate the angularly- averaged spec- 
trum and surface brightness profile. We have run standard test 
problems to test our 3D code and found that they are reproduced 
very well. 

Our modelling nevertheless has a number of significant 
simplifications. We post-process the output of a cosmological 
simulations. The simulations therefore do not take into account 
the dynamical effect of radiation pressure. Lya radiation pres- 
sure has been proposed to provide sufficient energy to launch 
an ou tflowing supershell for luminous sources (iDijkstra & Loebl 
120091) . However, the systems that we consider here do not have 
the necessary Lya luminosity for radiation pressure to be dy- 
namically significant. Our simulations also ignore the evolution 
of the density and velocity field over the light travel time of the 
Lya photons. We further only consider a very centrally concen- 
trated emissivity. For simplicity we inject photons at the centre 
of mass of the halo. The photons are created with line-centre 
frequency in the fluid frame of the gas. We will discuss these 
assumptions further in Section [STTI 



3 Lya RADIATIVE TRANSFER IN DLA HOST 
GALAXIES 

3.1 The Simulations of DLA host galaxies in a 
cosmological context by Tescari et al. 

We will make here us e of the cosmologica l hydrodynamic sim- 
ulations described in Tes cari et al ] d2009h . These simulations 
were aimed at reproducing the physical properties of the host 
galaxies of DLAs at z ~ 3. The numerical code is a modification 



of ga dget- 2, which is a parallel Tree-PM SPH code ( Springel 
120051) . The code is extensively described in iTornatore et al.1 
d2007l), where it was used to s imulate clusters of g alaxies; see 
also lTornatore et al] (120101) and lTescari et al.l (120101) . In addition 
to gravity and hydrodynamics, the other physical processes that 
are modelled in the simulation are: 

• Radiative cooling and heating, including a UV background 
produced by quasars and galaxies. 

• C hemical evolution, using the model of ITornatore et al.l 
(2007) which traces the following elements: H, He, C, O, Mg, S, 
Si and Fe. The contribution of metals is included in the cooling 
function. The release of metals from Type la and Type II super- 
novae, as well as low- and intermediate-mass stars, is followed. 

• Star formation, using a multiphase criterion. 

• An effective model for the ISM, whereby gas particles are 
treated as multiphase once they exceed a density threshold for 
the formation of cold clouds. These cold clouds are the precur- 
sors to stars. The neutral hydrogen fraction depends on the UV 
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background in the low density gas, and on the fraction f c of 
mass in cold clouds above the density threshold. 

• Galactic winds, using two different prescriptions. An 
energy-driven wind model is implemented in the form of a fixed 
velocity kick given to a chosen partic le, with a mass-loss rate 
proportional to the star formation rate. lTescari et al.l (12009b also 
investigated a momentum-driven wind model that mimics a sce- 
nario in which the radiation pressure of a starburst drives an 
outflow. In this model, the velocity kick scales with the velocity 
dispersion of the galaxy. These wind prescriptions are admit- 
tedly rather crude, and rely on phenomenological parameters 
that are poorly constrained either by observations or by more 
sophisticated modelling. This reflects the significant uncertainty 
regarding the influence of winds on galaxy formation. 



I Tescari et al investigated simulations with varying 

box size and numerical resolution for a range of different im- 
plementations of galactic winds. They chose the "SW" (Strong 
Wind) simulation with a box size of \0h~ l Mpc and 320 3 parti- 
cles, with a mass resolution of 3.5 x 10 5 h~ l M Q as their fiducial 
simulation. This simulation employed a strong (600 km s" 1 ), 
energy-driven wind and a Salpeter stella r IMF. Looking ahead 
to a favourable comparison with the lR08l emitters in a later sec- 
tion, we will instead focus more on the Momentum Driven Wind 
(MDW) simulation, which differs from the SW simulation only 
in the wind implementation. We will compare the results of both 
the MDW and SW simulation with two further simulations: a 
weak (energy driven, 100 km s -1 ) wind (WW) and a simula- 
tion w ith no wind (NW) which did not appear in ITescari et al.l 
d2009h . As we will see, the MDW and the SW simulations are 
somewhat similar, as are the WW and NW sim ulations. 

The simulations of T escari et al ] d2009h reproduce most 
observed properties of DLAs rather well. The observed inci- 
dence rate of DLAs is matched by the simulations, assuming 
that haloes below 10 9 /z _1 M o do not host DLAs. The observed 
column density distribution is also reproduced successfully for 
the SW and MDW simulations, while the total neutral gas mass 
in DLAs (^dla) is reproduced for MDW but underpredicted 
by a factor of about 1.7 at z = 3 in the SW simulation. The 
most problematic observable is the velocity width distribution 
of low-ionization associated metal absorption. The high veloc- 
ity tail of the distribution is significantly underpredicted in all 
simulations. The simulations do not produce enough absorp- 
tion systems with velocity width greater than 100 km s _1 . This 
is a com mon problem for numerical simulations of D LA host 
galaxies jPontzen et al.l 2008; iRazoumov et al.l l2008h . but see 
iHong et all d2010h and lcenl (I2010h which appear to be some- 
what more successful in this respect. We will come back to this 
point later. 



3.2 Three representative haloes 

We have selected three haloes at z = 3 from the simulations for 
detailed study. Some basic properties of the haloes are sum- 
marised in Table Q] The virial velocity and radius are calculated 
from the mass of the part icles grouped by the halo finder — see 
iMaller & Bullockl (120041) for the relevant formulae. Note that 
the properties of a given halo for the different wind model simu- 
lations are very similar. The fin al colum n l ists the Lya luminos- 
ity (L Lya ), calculated following IbH091 and lBHIOl as, 



where y c M 1 ^ 3 is the virial velocity of the halo. With this 
choice IBH09I and BH10 were able to reproduce the luminos- 
ity function of the lR08l emitters. The masses of the three haloes 
probe the range of masses predicted to contribute significantly 
to the incidence rate of DLAs / emitters per unit log mass 
(d 2 W/dX/dlog 10 M). The correspondin g lum inosities span the 
range of those observed in the survey of |R08|. 

The physical properties of these haloes (neutral hydrogen 
density, temperature and bulk velocity) were projected onto a 
regular 3D grid centred on the centre of mass of the halo. The 
cubes containing Haloes 1 and 2 have 128 3 cells, each with side 
length 3.125/T 1 comov.kpc, making the entire cube 400/T 1 co- 
mov.kpc across. The cube containing Halo 3 has 64 3 cells, each 
also 3.125/T 1 comov.kpc across, making the with of the entire 
cube 200/T 1 comov.kpc. 

We will first examine the simulation with the MDW imple- 
mentation of galactic winds, before comparing with the other 
simulations in later sections. Figure Q] shows the angularly- 
averaged radial mass flow rate, projected column density of neu- 
tral hydrogen and the Lya: surface brightness as a function of 
radius/impact parameter as well as the angularly- averaged spec- 
tral shape of the Lya- emission for the three halos. 

The rate at which Hi mass flows through a surface of con- 
stant radius r is shown in the upper left panel of Figure Q] nor- 
malised to a maximum \M\ = 1. Positive (negative) values cor- 
respond to outflow (inflow). The maximum M is 97, 8.5 and 1.5 
M©/yr, for Halo 1, 2 and 3 respectively. Note that, for the MDW 
implementation shown here, only in the most massive halo does 
the wind actually manage to globally reverse the gravitational 
infall. The mass flow shows a mixture of inflow and outflow, 
with inflow often dominating in the outer parts of the halo. 

The upper right panel of Figure Q] shows the average col- 
umn density for a sightline passing at a given impact param- 
eter from the centre of the halo. The thin dotted horizontal 
line indicates the minimum column density of a DLA, A^dla - 



10 



,20.3 



cm . As expected, the column density peaks near the 



VlOOkms" 1 / 



erg s 



-l 



(1) 



centre of the halo. The "steps" in the column density profiles 
of Haloes 2 and 3 show the influence of separate clumps of Hi 
at large radii. The neutral hydrogen column density exceeds the 
threshold for a DLA, which extends to a radius of 3.5, 1.5, and 
0.5 arcsec, respectively. 

The lower left panel shows the (observed) angularly- 
averaged spectrum, with each curve normalised to unity. The 
lower right panel shows the surface brightness profile. The lumi- 
nosity used to normalise each curve is given in the last column 
of Tabled] 

The surface brightness profiles peak at the centre of the 
halo, and generally follow the decline of the column density 
whilst being much smoother. The Lya- emission is scattered to 
similarly large radii as in the observed emitters. The haloes are 
surrounded by low surface brightness emission which extends 
to radii of several arcsec. The size o f the emitter, as measured 
to the surface brightness limit of the lR08l (shown by the hori- 
zontal dotted line), is generally larger than the cross-section for 
damped absorption. The more massive haloes are more spatially 
extended, while the smaller haloes display a more pronounced 
central peak. 

The angularly- averaged spectrum of the most massive halo 
is fairly symmetric with well separated blue and red peaks. The 
red peak is slightly stronger. The two low mass haloes have a 
more clearly dominant red and blue peak, respectively. Increas- 
ing central column density clearly leads to more diffusion in 
frequency space. We will see in the next section that the spec- 
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Halo Wind Mass Mass Hi Virial Velocity Virial Radius Lya Luminosity 

ID Model (M Q ) (Mo) ( kms" 1 ) (kpc) (arcsec) (xlO 40 erg s" 1 ) 



1 


momentum (MDW) 


7.70X10 11 


2.03xl0 10 


206.7 


77.5 


9.8 


882.7 




strong (SW) 


7.24X10 11 


1.38xl0 10 


202.5 


76.0 


9.6 


830.9 




weak (WW) 


8.15X10 11 


2.80xl0 10 


210.7 


79.0 


10.0 


935.3 




none (NW) 


8.12X10 11 


2.00xl0 10 


210.4 


78.9 


10.0 


931.5 


2 


MDW 


1.54X10 11 


3.36X10 9 


120.9 


45.3 


5.8 


176.5 




SW 


1.40xlO n 


1.07xl0 9 


117.2 


43.9 


5.6 


160.8 




WW 


1.70xlO n 


6.22xl0 9 


125.0 


46.9 


5.9 


195.1 




NW 


1.65X10 11 


4.76xl0 9 


123.7 


46.4 


5.9 


189.3 


3 


MDW 


1.46xl0 10 


2.74xl0 8 


55.1 


20.7 


2.6 


16.7 




SW 


1.36xl0 10 


1.30xl0 8 


53.8 


20.2 


2.6 


15.6 




WW 


1.58xl0 10 


8.10x10 s 


56.6 


21.2 


2.7 


18.2 




NW 


1.57X10 10 


9.12x10 s 


56.5 


21.2 


2.7 


18.1 



Table 1. Characteristic properties of three DM haloes of DLA host galaxies. 
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Figure 1. The angularly-averaged properties of the three DLA host galaxies described in Table 1, for the MDW wind model. The upper left panel 
shows dM/dt, which is the rate at which Hi mass is flowing through a surface of constant radius r, normalised to a maximum \M\ = 1. The maximum 
M (by which the curve are scaled) are 97, 8.5 and 1.5 M©/yr, for Halo 1, 2 and 3 respectively. The upper right panel shows the average column density 
for a sightline passing at a given impact parameter from the centre of the halo. The thin dotted horizontal line indicates the minimum column density 
of a DLA, jVdla = 10 203 cm -2 . The lower left panel shows the (observed) angularly-averaged spectrum, with the area under each curve normalised to 
unity. The lower right panel shows the surface brightness profile; the horizontal line shows the surface brightness limit of the R08 survey. The assumed 
total Lya luminosities (which normalise each curve) are given in the last column of Table[T] 



tral shape varies strongly between different directions and that 
there is a correlation between the dominance of inflow (outflow) 
and the domination of the blue (red) peak in the spectrum. We 
will leave a more detailed discussion until then. 



4 THE EFFECT OF GALACTIC WINDS ON Lya 
EMISSION 

We will now discuss the simulations of the intermediate mass 
halo (our fiducial Halo 2) in more detail with particular empha- 
sis on the effect of galactic winds on the physical properties of 
the gas in the halo and on the emission properties of the Lya 
emission. Halo 2 has a mass of about 1.5 x 10 10 M o and a virial 
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velocity of about 120 km s" 1 . Its p roperties are similar to the 
brightest of the emitters in the [R08 survey, and its mass coin- 
cides with the peak of our preferre d mod el for the DLA mass 
function (d 2 W/dX/dlog 10 M) from lBHlol and thus should rep- 
resent a typical DLA host galaxy. 



4.1 Angularly- Averaged Properties 

Angularly- averaged properties of simulation based on Halo 2 
for each of the wind models are shown in Figure [2] The legend 
in the top right panel shows the line- styles of the different wind 
models 

The top left panel shows the Hi number density as function 
of radius. Note that the average cosmic number density of hy- 
drogen at z ~ 3 is 1.2 x 10" 5 cm -3 , while the volume- averaged 
neutral fraction of hydrogen in the IGM is < 10" 5 . A small de- 
crease in the number density at small radii is the result of the 
centre of mass of the neutral hydrogen being offset from the 
centre of mass of the halo. The number density profile shows 
the effect of several clumps of Hi. The most dramatic difference 
between the models is in the centre of the haloes, where the 
presence of a strong wind can reduce the gas density substan- 
tially. 

The top right panel shows the averaged column density for 
a sightline passing at a given impact parameter from the cen- 
tre of the halo. The thin dotted horizontal line indicates the 
minimum column density of a DLA. The column density is a 
decreasing function of radius, with sightlines containing DLAs 
probing the centre of the halo. The "clumps" noticeable in the 
number density distribution, cause steps in the angle- averaged 
column density distribution and contribute to the Lyman-Limit 
System cross-section. 

The next panel (row 2, column 1) shows the (Hi) density- 
weighted radial velocity profile, 



vm(r) = (Vb ■ r) = 



J Q 71 £ (v h ■ r) p H i(r, 6>, 0) sin<9 d<9 d0 
jf £ PHi(r,0,0)sin0d0d0 



(2) 



The plot reveals the presence of a mixture of inflow (v r negative) 
and outflow, with the innermost parts of the halo tending to show 
outflow for the stronger winds (SW and MDW), and inflow for 
the weaker winds (WW and NW). 

The typical velocities of inflow are 50 - 100 km s -1 , with 
the largest velocities being comparable to the virial velocity. 
There is a general trend toward lower velocities in the centre of 
the halo. The velocity profile is, however, far from smooth. The 
total velocity of the gas (v = V(N>I 2 )X by comparison, varies 
between 100-150 km s -1 for the stronger winds, and 140 - 200 
km s" 1 for the weaker winds, with a very slight preference for 
larger velocities at the centre of the halo. 

The panel to the right shows the volume-weighted ra- 
dial velocity profile (vy). For the diffusion of Lya photons, 
volume-weighted velocities are actually more relevant than 
density/mass-weighted velocities as the photons tend to find 
low-density paths of escape, and so the velocity in the under- 
dense regions is as important or more important as that in high- 
density regions. We see a clear difference in vy between the 
stronger and weaker wind models. The strong wind models have 
outflows that dominate the kinematics of the gas along most of 
the escape paths of the photons. For weak wind models, inflow 
dominates. The next panel (row 3, column 1) shows dM/dt for 
Hi, which is related to the density-weighted radial velocity as 



dM/dt = 47Tr 2 vu l (r)p^ i (r), where pj^ is the average neutral hy- 
drogen density at radius r. Notice that the flows in the innermost 
parts of the halo, while being of modest velocity, carry the most 
significant amounts of mass. 

Comparison of vn^r), vy and dM/dt is quite instructive. 
The relative smoothness of vy indicates that a small number of 
dense clumps dominate the mass flow, moving through a rela- 
tively smooth background whose kinematics are shaped by the 
winds. For example, the MDW model (solid black) shows a 
~ 100 km s" 1 inflowing clump at 2 arcsec that carries a sig- 
nificant amount of mass. However, vy is still positive, showing 
that the wind is filling the remaining space and outflowing re- 
gardless. Lya photons escaping along low column density paths 
will generally not "see" such dense clumps. 

The panel to the right shows the rotation velocity of the 
gas, averaged over cylindrical shells, which is approximated as 
follows. If the gas in the halo were in a rotating disk with a rota- 
tion velocity depending only on the cylindrical radius, then the 
angular momentum within a cylindrical shell R,R + dR would 
be, 



ILrotl = v TOt (R) R^rrii 



(3) 



where the sum is over the cells inside the shell. Thus for each 
R, we divide the angular momentum by R times the mass in the 
shell. As with um(r), this distribution is far from smooth due to 
the influence of dense clumps of gas. The velocity drops from 



~ 150 km s to ~ 50 km s as we move outwards. The ve- 
locity is generally higher in the weaker wind models, possibly 
indicating that the stronger winds disrupt or delay the settling of 
gas into a disk. We can also calculate the fraction of the kinetic 
energy at each radius which is in the form of rotational kinetic 
energy. This quantity peaks at ~ 0.7 at the centre of the halo, 
decreasing to ~ 0.2 in the outer parts. This shows that the ap- 
proximations we made in calculating v YOt are most accurate at 
the centre of the halo, as expected. 

The bottom left panel shows the angularly- averaged spec- 
trum, of the Lya emission. The thick black curve shows the 
spectrum as a function of the rest-frame velocity offset (i> p h = 
cAA/Aq). All curves are normalised to have unit area. The 
weaker wind models show a c lassic dou ble-peaked profile 
faarringtonl 1 19731 : iNeufeldl 1 19901 : see also lUrbaniak & Wolfd 
|_198JJ), while the stronger wind models show a more prominent 
red peak, indicating the influence of the outflowing gas. 

We also see that the stronger galactic winds — somewhat 
counterintuitively — result in a narrower spectral distribution. 
This indicates that it is the high central column density, more 
than the higher outer gas velocities, that are the dominant factor 
in establishing the typical escape frequency of Lya photons. 

The thick black curve in the bottom right panel shows the 
angularly- averaged surface brightness profile. The luminosity 
that normalises each curve is given in the last column of Tablep] 
The vertical line indicates the virial radius of the halo, whil e the 
horizontal line shows the surface brightness limit of lR08l . The 
surface brightness profile shows a central peak with a tail that 
flattens out as the column density does the same. The emission 
is more extended in the stronger wind models. There are two 
reasons for this. The lower central column density means that 
photons are in general closer to line centre when they encounter 
the outer parts of the halo, where higher wind velocities can shift 
the photons sufficiently back towards line centre to be scattered 
by rather low column density gas. Comparing the radial profile 
of Hi column density and Lya surface brightness for the differ- 
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Figure 2. Angularly-averaged properties of Halo 2, for the wind prescriptions as shown in the legend (MDW: momentum driven wind, SW: strong 
wind, WW: weak wind, NW: no wind). The velocities and dM/dt are all for Hi. Vertical dashed lines indicate the virial radius. For a detailed discussion, 
see Section |4~T1 
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ent wind models suggest that this is more important than the the 
Hi column density in the outer parts of the halo. 

We have not plotted the (Hi weighted) temperature profile, 
which for all models decreases gently with radius from around 
10 5 K to 10 4,3 K. The central temperature is slightly higher in 
the weak wind models. 

We can now c ompare the properties of the haloes to the 
modelling of lBHld Unlike in our previous analytical modelling 
the column density in our simulations here has a pronounced 
core at very small radii. At larger radii the radial profile of the 
Hi column density in the MDW model appears rather similar 
even though much larger samples of simulated haloes would 
be needed t o make a more mea ningful comparison; we refer 
the reader to lTescari et al.l j2009h for a quantitative comparison 
with the observed DLA column density distribution. The radial 
velocity profiles in our full cosmological simulations are obvi- 
ously much more complex than the simple power law profiles 
we had assumed in our previous work, although there is a slight 
trend toward lower velocities at the centre of the halo. There 
was obviously also no possibility to account for the 3D effects 
of any clumping in our previous ID modelling. 

4.2 2D Views of the Halo 

Figure [3] shows maps of the Hi column density distribution and 
Lya surface brightness and 2D spectra for Halo 2 for the MDW 
model looking along the x-axis from -oo (left) and +00 (right). 
Note that some of the axes have been flipped so that the left and 
right columns have their axes pointing in the same direction. 
The dashed black curve shows the virial radius. 

The top panels show the Hi column density along a line 
of sight parallel to the x-axis. Each panel corresponds to a line 
of sight that passes through the nearer half of the box. In other 
words, the colours in the legend to the right of each Figure rep- 
resent the column density of neutral hydrogen that the photon 
(emitted at the centre of the halo) would have to pass through 
if it were to reach the observer directly without any scattering. 
The solid line in the top two panels is a contour representing 
A^dla - 10 20 ' 3 cm -2 , calculated for sightlines passing through 
the entire box; it is thus the same in the left and right top panels. 
A sightline passing through the region inside the contour would 
encounter a DLA. 

The column density distribution reflects the irregular dis- 
tribution of gas, with a mixture of thin filaments and relatively 
isolated clumps of neutral hydrogen. The differences between 
the panels show the asymmetry in the distribution of gas. The 
cross-section for damped absorption is composed of a large cen- 
tral clump accompanied by a number of isolated clumps, scat- 
tered around the halo. The column density of individual clumps 
of neutral hydrogen drops rapidly below ~ 10 17 cm -2 , due to 
the effect of the ionising UV background on particles below the 
'cold cloud' threshold. 

The next two panels show an image of the halo that is 
coloured according to the surface brightness of the Lya emis- 
sion, calculatecQ assuming the luminosity as described pre- 
viously and given in Table Q] The solid contour represents 
the surface brightness limit of the lR08l survey, So = 10" 19 
erg/s/cm 2 /arcsec 2 . 

The surface brightness maps show that the Lya emission 
generally traces but is less clumpy than the underlying neutral 

1 See equation (20) of lLaursen, Razoumov, & Sommer-La rsen (2009). 



hydrogen distribution. In particular, the So contour encloses a 
much larger area than the DLA contour both by being larger 
and also by having fewer "holes". The Lya photons are scattered 
effectively out to large radii and lower column densities. The So 
contour in the surface brightness is roughly similar to that of the 
Nhi ~ 10 18 cm -2 contouinin the Hi column density plot. 

Our central Lya source is able to illuminate clouds of neu- 
tral gas that are detached from the central clump. The Lya pho- 
tons tend to find low-density paths of escape. For example, the 
photons face a column density of Hi which is on average ~ 7 
times larger in the -z direction than in the +z direction, and 
the corresponding observed Lya flux is on average ~ 3 times 
smaller. For the x direction, as shown in the Figure, the differ- 
ence in the observed flux between two opposite viewpoints is 
around 70%. 

The bottom four panels are 2D spectra, calculated us- 
ing a "slit" placed along either the y- or z-direction (labelled 
in the vertical ax is) and with a wid th of 2 arcsec, the same 
width as used in iRauch et all d2Q08h . The contour shows the 
spectral intensity limit of the Ir08| survey, S ^ « 4 x 10" 20 
erg/s/cm 2 /arcsec 2 /A. 

The 2D spectra show a dominant blue peak in the -x di- 
rection, and the opposite trend when viewed from +x; see the 
discussion of Figure Hlbelow. As expected, the separation of the 
peaks is largest where the column density is largest. The dom- 
inance of one peak over the other tends to increase as we look 
away from the centre of the halo. This is most likely because the 
probability of a photon scattering in the outer parts of the halo 
is quite sensitive to the velocity of the gas. 

We illustrate the dependence of the emergent spectrum 
and surface brightness profile on the viewing angle in Figure 
|H which also shows the angularly- averaged value (thick black 
line). The thin coloured lines show the emergent spectrum (cal- 
culated by collapsing the slit along its spatial direction, or by 
averaging over the 2D image) for 6 different viewpoints along 
the principal axes. Remember that the observer sees only one 
of these spectra - the angularly- averaged quantities are not ob- 
servable. The spectra show significant variations, both in which 
peak dominates and in total amount of light observed. The sep- 
aration of the peaks shows relatively little variation and the sur- 
face brightness profiles are very similar, especially in the outer 
parts of the halo. 



4.3 Haloes 1 and 3 

The angularly- averaged properties of Halo 1 are shown in Fig- 
ure[5] The panels show the volume- averaged radial velocity pro- 
file, projected column density of neutral hydrogen, the emer- 
gent Lya spectrum and the Lya surface brightness as a function 
of radius/impact parameter. The four different wind models are 
shown on each plot. 

The radial velocity profile shows the effect of the outflow- 
ing winds in the low density gas. The winds provide escape 
paths filled with gas that is outflowing at ~ 100 - 150 km s" 1 , 



2 Varying between ~ 10 17 — 10 19 cm -2 . Note that this is not the 
same as the mean integrated column density along the final escape 
flight of the photon, which we have also calculated explicitly us- 
ing our code and which can be approximated analytically by N GSC « 
10 20 cm -2 (^peak/500 km s -1 ) 2 , where i; pea k is the peak of the emission 
spectrum. 
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Figure 3. 2D images and spectra for Halo 2 in the MDW simulation. The left (right) column is as viewed from x = -00 (+00). The upper two panels show 
the Hi column density for the closer half of the simulation box. The dashed circle shows the virial radius. The contour represents A^dla = 10 20 3 cm -2 , 
for a sightli ne pa ssing through the entire box. The next two panels show the observed Lya surface brightness. The contour shows the surface brightness 
limit of the R08 survey, So = 10~ 19 erg/s/cm 2 /arcsec 2 . The lower four panels show 2D spectra, c alculated using a "s lit" placed along either the y- or 
z-direction (labelled in the vertical axis) and with a width of 2 arcsec, which is the same width as in lRauch et al.l d2008h . The contour shows the spectral 
intensity limit of the lR08l survey (Rauch, M., private communication), which is S ^ ~ 4x 10~ 20 erg/s/cm 2 /arcsec 2 /A. The Figure is discussed further in 
Section l4~2l 
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Figure 4. The emergent spectrum (left) and surface brightness profile (right) for Halo 2, showing the effect of the viewing angle. The thick black line 
is the angularly-averaged spectrum, which is the same as in Figure [T] The thin coloured lines show the emergent spectrum (calculated by collapsing the 
slit along its spatial direction, or by averaging over the 2D image) for the 6 different viewpoints as shown in the legend. The Figure is discussed further 
in Section l4~2l 
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Figure 5. Angularly-averaged properties of Halo 1, for the wind prescriptions shown in the legend (upper left), as discussed in Section 1431 The upper 
left panel shows the volume-averaged radial velocity profile. The upper right panel shows the Hi column density as a function of impact parameter. 
The lower left panel shows the spectra, normalised to unit area. The lower right panel shows the surface brightness distribution as a function of impact 
parameter. 



while in the weaker wind cases the average velocity is closer to 
zerc0. 

The column density profile of Halo 1 is somewhat 



smoother than that of Halo 2 and 3, and as expected the DLA 
cross section is significantly larger. Once again we see that the 
winds have lowered the central column density of the halo. 



3 The fact that the volume-averaged radial velocity is slightly positive 
in the no-wind simulation is most likely due to the centre of mass of 
the baryons being offset from the centre of mass of the halo i.e. the 



baryons have overshot the bottom of the potential well. Heating of the 
hot, rarefied ISM by the UVB may also contribute. 
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For Halo 1 the spectra show again the classic double- 
peaked profile with a general trend toward domination of the 
red peak for haloes with stronger winds. The surface brightness 
profiles are rather similar: they peak at the centre, and decrease 
smoothly as the column density drops. The weaker wind mod- 
els are as expected more centrally peaked than the stronger wind 
models. 

The angularly- averaged properties of Halo 3 are shown in 
Figure(6] The gas in the stronger wind models is more extended. 
Volume-weighted all of the models are dominated by infall at 
small radii, with only the SW model showing a very mild out- 
flow at large radii. The MDW wind clearly decreases the average 
infall rate, but not to the extent that the infall is reversed. 

The angularly- averaged spectrum of Halo 3 is a particu- 
larly nice illustration of inflowing gas at the centre of the halo 
producing a more prominent blue peak in the spectrum. In spite 
of this, there are viewing angles for which the spectrum is dom- 
inated by a red peak — see Figure [9] below. The surface density 
profile reflects the column density profile — the inner, dense 
core results in a central peak whose rapid decline gives way to 
a shallower plateau as we reach the outer, diffuse background. 

Figures [7] [8] and [9] show a comparison of the 2D maps and 
spectra for the MDW (left) and WW (right) wind prescriptions 
for the three haloes. Note that the left panel of Figure [8] is the 
same as the left panel of Figure [3] and is reproduced in this 
Figure for convenience. As noted above, the emission is more 
extended in the stronger wind models, and stronger winds pro- 
duce a narrower spectral distribution (note the difference in the 
horizontal scale of the 2D spectra). 



5 DISCUSSION 

The dependence of the spectra, in both the line shape and the 
total observed surface brightness, on the viewing angle is very 
striking and is obviously very important for the interpretation of 
observed (spatially extended low surface brightness) Lya emit- 
ters. It suggests, g.ff. , that the large variety of line shapes in the 
observed sample of lR08l is at least partially an orientation ef- 
fect. The dependence of surface brightness on orientati on fur- 
ther suggest that the "duty cycle" that we introduced in iBHIOl 
may also be interpreted as an orientation effect, with only a cer- 
tain fraction of haloes seen at an orientation where the observed 
flux exceeds the detection limit. A rather large sample of haloes 
would be required to determine the probability distribution of 
the observed Lya flux from a halo of given mass and intrinsic 
luminosity. As briefly mentioned but not investigated in detail 
here, there should, however, obviously be also a variation of the 
spectral shape and surface brightness distribution with time in 
particular when a galaxy undergoes a major merger and/or star- 
burst. 



5.1 Limitations of the modelling 

We will discuss now again in more detail some of the limita- 
tions of our modelli ng. Firstly, the failure of the simulations of 
iTescari et al ] d2QQ9h to reproduce the high velocity tail of the 
velocity width distribution of the associated metal absorption in 
DLAs suggests that higher outflow velocities, a more extended 
spatial distribution of the neutral gas, and/or a less clumpy and 
more spatially extended distribution of low-ionisation metals is 
required from our galactic wind implementations. The effect of 



higher outflow velocities on the Lya emission depends on ex- 
actly how the velocity would be increased. If the bulk velocity 
in the centre of the halo is increased, then this would accelerate 
frequency- space diffusion resulting in a more centrally peaked 
surface brightness profile. On the other hand, if the bulk veloc- 
ity is increased primarily in the outer parts of the halo, then this 
would have a minimal effect, perhaps shifting photons back to- 
ward line-centre and thereby increasing the apparent size of the 
emitter. 

We have also not attempted to model the effects of dust on 
the Lyar radiation. DLAs are known to be among the most metal- 
poor systems in the high r edshift universe, and are t hus expected 
to be relatively dust free dPontzen & Pettinill2009h. apart from 
perhaps the most metal-rich systems dFvnbo et al . 2010a,rJ. The 
long random walk of a typical Lya photon makes them particu- 
larly susceptible to even a small amount of dust. 

Another worry could be our use of a rather coarse regular 
grid. Small-scale inhomogeneities in the Hi density are, how- 
ever, not likely to have much effect in the absence of dust. As 
we have seen the Lya emission profile smoothes them over. At 
most, an irregular edge of an Hi region may lead to a shallower 
surface brightness profile. Inhomogeneities in the bulk velocity 
would have a somewhat greater effect, increasing the frequency 
diffusion. These inhomogeneities could be modelled as a turbu- 
lent velocity, which has a similarly small effect as raising the 
temperature. 

The fact that we have simply injected Lya photons at the 
centre of the halo i s a rather strong albei t observationally mo- 
tivated assumption. IWolfe & Chenl ([2006) placed limits on the 
spatial extent of in situ star formation in DLAs, concluding that 
ext ended so urces (r adius > 3kpc, 0.4 arcsec) were very rarfl 
As lR08l and lBH09l have argued, Lya emission in the systems 
we are considering is most likely to be powered by star for- 
mation, suggesting that their large observed sizes are the re- 
sult of radiative transfer through the surrounding neutral gas. 
Our main focus here has been on whether such a scenario is in- 
deed plausible. The fact that the Hi density peaks very near the 
centre of the halo means that the majority of photons should 
be emitted near the centre of the halo, as we have assumed. 
Making the intrinsic Lya emission more extended would ob- 
viously make the observed emission also more extended. A 
more realistic Lyor emissivity would need to consider the in- 
teraction of UV radiation from star formation with the sur- 
rounding Hi, as well as the possible contr ibution from cool- 
ing radiation jFaucher-Giguere et "all l2010h and fluorescence 
due to the UV background and qua sars dCantalupo et al1l2005l ; 
ICantalupo. Porciani. & Lilly! l2008h . However, fluorescence is 
expec ted to contribute little to the Lyor emission from DLAs 
jR08h , and cooling is likely to only be significant in massive 
haloes, which aren't expected to host the bulk of the DLA pop- 
ulation. While cooling radiation may not account for the bulk 
of DLA lya luminosity, it could still be a relevant factor in 
determining the observed size of the emitters, as the surface 
brightness of the cooling radiation may drop less rapidly to- 
wards large radii than the radiation scattered outwards from 
central star-forming regions toiikstra. Haiman. & _Spaansll2006l : 
iDiikstra & Loeb|[2009l : lFaucher-Giguere et alj|2010h . 



4 Some extended star - format ion has recently been detected 
iRafelsk i. Wolfe, & Chenl 120101) — as discussed previously, we 
will investigate the effect of spatially extended Lya sources in future 
work. 
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Figure 6. Angularly-averaged properties of Halo 3, for the wind prescriptions shown in the legend (upper left), as discussed in Section POl The panels 
plot the same quantities as in Figure [5] 



We have considered here a rather small sample of haloes, 
for several reasons. Firstly, it allows us to discuss some of the 
many factors relevant for the complicated interpr etation of the 
Lya e mission at length. Also, the simulations of iTescari et al.1 
(120091) used a box size that was too small to produce large num- 
bers of haloes with M v > 1O 1O ' 5 M . This was necessary to give 
an acceptable resolution. Our w ork her e should thus be seen 
as complementary to the work of lBHld Many of the assump- 
tions made there are borne out in the simulations described here. 
The Lya emitters have considerable spatial extent, and are larger 
than the corresponding DLA cross-section as Lya photons can 
be effectively scattered by gas with column densities less than 
that of a DLA; the spectra are often double-peaked, with in- 
flow/outflow producing a dominant blue/red peak; the radial ve- 
locity and mass flow rate are quite moderate toward the centre 
of the halo, increasing the spatial extent of the Lya emission. 

5.2 Relation to LBGs 

Observationally well studied LBGs at similar redshift - presum- 
ably the more luminous and more massive cousins of typical 
DLA host galaxies - offer the opportunity to investigate their 
Lya emission in considerably more detail, though a compari- 
son may be hampered due to the increased role of dust and the 
fact that galactic winds are li kely to be more pow erful in these 
rapidly star-forming galaxies. Fsteidel et aTI(l201Qh have summa- 
rized the properties of the observed Lya emission in LBGs. 
Lya is seen in emission and absorption with a wide range of 
EW. When observed in emission, Lya typically shows a dom- 



inant red peak with offsets in the range of 200-800 km s 
with respect to the systemic redshift as defined by the nebu- 
lar emission lines. The medium offset of the Lya emission is 
485 + 175 km s" 1 with a typical width of 500 km s" 1 (e.g. 
IOuideretal.ll2009h . While the Lya emission is offset towards 
the red, the interstellar absorption lines are offset to the blue 
with a mean of about -250 km s" 1 and a tail of absorption of- 
ten extending out to -800 km s . In particular, the blueshift 
of the interstellar absorption line has been convincingly inter- 
preted as being due to outflowing gas moving towards the ob- 
server ( Verham me. Schaerer. &MasellJl2006h . The most mas- 
sive of our three haloes has a mass which is approaching that 
inferred for LBGs DM haloes so it is interesting to see how our 
simulations compare. It is the momentum wind driven model 
with its pronounced red peak which appears to have the most 
similar spectral distribution. The peak of the angular averaged 
emission in the MDW model is at about +400 km s -1 . The sec- 
ondary blue peak is, however, co nsiderably stronger t han in the 
composite spectrum presented in ISteidel et al.l d201Qh . perhaps 
indicating that the central column density is somewhat too large 
and that real galactic winds in LBGs are more efficient in reduc- 
ing the central optical depth. In fact, most LBGs do not show 
an intrinsic damped absorption system in their spectra and the 
outflow velocities indicated by the interstellar absorption lines 
appear to be higher than that of the outflowing gas in our simu- 
lations. 
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Figure 7. 2D images and spectra for the MDW (left) and WW (right) 
in Figure[3] This Figure is discussed further in Section 1431 
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prescriptions for Halo 1 viewed from x = -oo. The panels are the same as 
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Figure 8. 2D images and spectra for the MDW (left) and WW (right) 
in Figure[3] This Figure is discussed further in Section |4~3l 
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Figure 9. 2D images and spectra for the MDW (left) and WW (right) wind prescriptions for Halo 3 viewed from x = - 
in Figure[3] This Figure is discussed further in Section |4~3l 
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5.3 A more detailed comparison with the Rauch et al. 
emitters 



We come now to a more detailed compa rison to the R08| emit- 
ters. As discussed in lR08l (see also lBHlOk for 12 of the 27 spec- 
tra only a single emission peak is visible while six/three of the 
spectra show a weak secondary blue/red counter-peak. The re- 
maining spectra are extended in frequency space without a clear 
peak structure. The widths of the spectral peaks ranges from 
~ 250- 1000 km s -1 . The surface brightness profiles are pre- 
dominantly centrally peaked with wings that often extend well 
beyond the Gaussian core of the PSF. This is particularly true 
of the brightest sources, while the fainter sources are more dif- 
ficult to characterise due to the noise. As discussed in detail in 
section |4l for the weak and no wind model the spectral distribu- 
tion shows rather symmetric double-peaked profiles. The lR08l 
emitters appear instead to be similar to the two models with the 
stronger winds. Of the two models with strong winds, the MDW 
model shows the more pronounced preference for a dominant 
red peak and generally the contrast between the two peaks is 
largest in this model. This model is thus probably the one which 
comes closest to the observed emitters. As discussed, the main 
difference between the models with the different wind imple- 
mentations is due to the reduction of the central column density 
of neutral hydrogen due to the effect of the galactic wind. Note 
that most of the R08 emitters are only just above the detection 
threshold and offer thus a somewhat limited dynamic range. A 
weaker second peak would probably not be detected in many 
of the observed single-peak emi tters. There is, however, a no- 
ticeable difference: some of the lR08l emitters have more than 
one peak in the spatial direction suggesting a more spatially ex- 
tended source of ionizing photons perhaps with several peaks. 
As discussed previously, in our simulation the photons were all 
injected at the centre of the DM halo. As is apparent from the Hi 
column density maps, the Hi distribution is actually also rather 
clumpy. In a more quantitative study of a larger sample of haloes 
it would thus certainly be worthwhile to test if the correspon- 
dence between observation and simulations can be improved by 
injecting Lya photons in several Hi density peaks. 

5.4 Implication for the detectability of faint high-redshift 
Lya emitters 

The velocity shift of Lya emission due to resonant scat- 
tering is very important for the detectabili ty of Lya emit- 
ters at high redshift (e.g. ISantosI 120041) . As hydrogen 
becomes increasingly neutral at high redshift, the opti- 
cal depth to Lya scattering due to the neutral IGM be- 
tween the emitters and the observers steadily increases. 
As di scus sed e.g. by Zheng et al. (2009), Diikstr a~& Wvithd 
(l2010h and IE aursen, Sommer-Larsen, & Razoumov (2010), the 
effective transmission thereby depends sensitively on the op- 
tical depth to Lya scattering and the details of the gas ve- 
locities within the emitter and the intervening IGM. Ve- 
locity shifts of about 500 km s" 1 as observed in z ~ 3 
LBG/LAEs would strongly increase the tra nsmission which 
for appreciable shifts to the red scales as ( Miralda-Escude 
ll998l : lDiikstra. Lidz. & Wvithd2007tlmikstra & Wvithel2010h . 
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much fainter IR08I emitters still suggests shifts of about 300 
km s -1 which is good news for the detectability of faint Lya 
emitters at high redshift perhaps even before reionization is 
complete. Note, however that in our model the shift towards the 



red decreases rather than increases with increasing strength of 
galactic winds. As discussed in Section[4] this reduced shift to- 
wards the red is due to the decreased central column density of 
neutral hydrogen in the models with the stronger galactic winds. 
For a good estimate of the expected velocity shift it should thus 
be more important to predict the correct column density than the 
exact outflow velocity of the gas. 



6 CONCLUSIONS 

We have used here realistic numerical simulations of DLA host 
galaxies including the effect of galactic winds to investigate the 
spatial and spectral distribution of low-surface brightness Lya 
emission assuming that the emission is powered by star forma- 
tion at the centre of the corresponding DM haloes. Our main 
results can be summarised as follows. 

• The haloes contain a mixture of inflowing and outflow- 
ing gas. As a result, the angularly- averaged spectrum typically 
shows two peaks, with the relative strength of the red (blue) 
peak being a reflection of the relative contribution of outflow 
(inflow) as characterized by the volume- weighted radial veloc- 
ity of Hi. The separation of the two peaks is mainly governed by 
the central Hi column density. Line-centre photons can escape 
due to inhomogeneities in the gas density providing low-density 
paths of escape and bulk velocities allowing rapid diffusion in 
the fluid frame frequency space. The different wind implemen- 
tations lead to significantly different central Hi column densities 
as well as different radial motions of the gas, which are reflected 
in different spectral shapes of the Lya emission. 

• A comparison of the 2D column density and surface bright- 
ness images shows that the Lya emission region is larger and 
smoother than the cross-section for damped absorption. Lya 
photons escaping at large radii illuminate regions of protogalax- 
ies that would be probed by absorption line spectra with col- 
umn densities down to A^hi ~ 10 18 cm -2 and below. The cen- 
tral source can effectively light up outlying clumps of neutral 
hydrogen. Asymmetries in the Hi distribution can dramatically 
affect the observed surface brightness, with dense clumps cast- 
ing shadows as photons are reflected off their surface. A central 
source of photons can illuminate diffuse, highly ionised gas in- 
side and beyond the virial radius at very faint surface brightness 
levels. 

• The 2D spectra show considerable variety for the same halo 
viewed from different angles. A typical line profile is double- 
peaked with one of the peaks dominating. The separation of 
the peaks decreases with increasing distance from the central, 
dense regions of the halo. The dominance of the red or blue 
peak changes with viewing angle, and is affected by the veloc- 
ity of the gas relative to the observer. The maximum spectral 
intensity is approximately constant across the different haloes at 
S a « 10" 17-5 erg/s/cm 2 /arcsec 2 /A, with the decrease in luminos- 
ity with mass being countered by reduced spatial and frequency 
diffusion. 

• The angularly- averaged surface brightness profile is most 
sensitive to the column density of the gas at the centre of the 
haloes which in turn is very sensitive to the feedback from galac- 
tic winds. The more efficient feedback implementations result 
in reduced column densities at the centre and therefore reduced 
diffusion in frequency space and narrower spectral profiles. The 
maximum blue/redshift in our emitters is also primarily sensi- 
tive to the central HI column density and not the velocity of 
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the wind. For objects with column densities and wind veloc- 
ities similar to those in our simulations, galactic winds may 
in fact reduce rather than increase the visibility of Lya emit- 
ters during the epoch of reionization when the hydrogen in the 
IGM is not yet (fully) ionized and Lya emission from galaxies 
is strongly suppressed due to the damping wing of the IGM's 
Gunn-Peterson trough. 

• For the momentum-driven wind model (and to a some- 
what lesser extent also the strong wind implementation), our 
simulated em itters show encouraging agreement with the prop- 
erties of the Ir08| emitters and thus further corroborate the 
suggestion that these are the long- sought host population of 
DLAs. This preference for the momentum-driven wind mod- 
els also appears to result in better agreement with a range of 
other observed properties of h igh-redshift-galaxies and the IGM 
dOppenheimer & Pavel l20Q8h . There are, however, at the same 
time differences between observed and simulated emitters. The 
simulated spectra do not show quite the preference for a dom- 
inant red peak the observed emitters appear to suggest. This is 
partially due to the fact that we have not modeled the impact of 
the intervening Intergalactic Medium which will reduce the ob- 
served emission in the blue peak relative to the red peak. Our in- 
vestigation nevertheless appears to suggest that somewhat larger 
outward motion are required than present in our MDW simula- 
tions where especially in the smaller haloes the feedback due 
to the galactic winds is not quite sufficient to globally reverse 
the infall of gas. This may well be related to the fact that the 
simulations also struggle to reproduce the tail of high velocity 
widths for the associated low-ionization absorption of DLAs. 
The existence of several clumps in the observed Lya emission 
further suggests that the emissivity of ionizing photons is more 
spatially extended than the central injection we have assumed in 
our simulations. 

The study of the effect of galactic winds in a cosmologi- 
cal galaxy formation simulation presented here contributes fur- 
ther to unravelling the complexities of the Lya emission from 
high redshift galaxies. The rather strong dependence of the Lya 
emission on the details of the wind implementation should al- 
low spatially extended low-surface brightness Lyor emission to 
be turned into an important diagnostic tool for studying the in- 
ternal dynamics of galactic winds as well as their effect on the 
interstellar and circumgalactic medium of high-redshift galax- 
ies. 
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